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Abstract. This study is motivated by the unusually low 
number of magnetic clouds (MCs) that are strictly identified 
within interplanetary coronal mass ejections (ICMEs), as ob- 
served at 1 AU; this is usually estimated to be around 30% or 
lower. But a looser definition of MCs may significantly in- 
crease this percentage. Another motivation is the unexpected 
shape of the occurrence distribution of the observers’ “clos- 
est approach distances” (measured from a MC’s axis, and 
called CA) which drops off somewhat rapidly as |CA| (in % 
of MC radius) approaches 100%, based on earlier studies. 
We suggest, for various geometrical and physical reasons, 
that the |CA|-distribution should be somewhere between a 
uniform one and the one actually observed, and therefore the 
30% estimate should be higher. So we ask, When there is a 
failure to identify a MC within an ICME, is it occasionally 
due to a large |CA| passage, making MC identification more 
difficult, i.e., is it due to an event selection effect? In attempt- 
ing to answer this question we examine WIND data to obtain 
an accurate distribution of the number of MCs vs. |CA| dis- 
tance, whether the event is ICME-related or not, where ini- 
tially a large number of cases (N= 98) are considered. This 
gives a frequence distribution that is far from uniform, con- 
firming earlier studies. This along with the fact that there 
are many ICME identification-parameters that do not depend 
on |CA| suggest that, indeed an MC event selection effect 
may explain at least part of the low ratio of (No. MCs)/(No. 
ICMEs). We also show that there is an acceptable geomet- 
rical and physical consistency in the relationships for both 
average “normalized” magnetic field intensity change and 
field direction change vs. |CA| within a MC, suggesting that 
our estimates of |CA|, Bo (magnetic field intensity on the 
axis), and choice of a proper “cloud coordinate” system (all 
needed in the analysis) are acceptably accurate. Therefore, 
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the MC fitting model (Lepping et al., 1990) is adequate, on 
average, for our analysis. However, this selection effect is not 
likely to completely answer our original question, on the un- 
expected ratio of MCs to ICMEs, so we must look for other 
factors, such as peculiarities of CME birth conditions. As a 
by-product of this analysis, we determine that the first order 
structural effects within a MC due to its interaction with the 
solar wind, plus the MC’s usual expansion at 1 AU (i.e., the 
non-force free components of the MC’s field) are, on average, 
weakly dependent on radial distance from the MC’s axis; that 
is, in the outer reaches of a typical MC the non-force free ef- 
fects show up, but even there they are rather weak. Finally, 
we show that it is not likely that a MC’s size distribution sta- 
tistically controls the occurrence distribution of the estimated 
| C A | s . 

Keywords. Interplanetary physics (Interplanetary magnetic 
fields; Solar wind plasma) 


1 Introduction 

A magnetic cloud (MC) was originally defined empirically in 
terms of in-situ spacecraft measurements of magnetic fields 
and particles in the interplanetary medium, viz., it is a re- 
gion in the solar wind having: (1) an enhanced magnetic field 
strength, (2) a smooth change in field direction as observed 
by a spacecraft passing through the MC, and (3) a low pro- 
ton temperature (and low proton plasma beta) compared to 
the ambient proton temperature (Burlaga et al., 1981, 1990; 
Klein and Burlaga, 1982; Burlaga, 1995). MCs are also un- 
derstood tacitly to be large structures averaging about 20 h 
in duration at 1 AU for the better examples (e.g., Lepping et 
al., 2006). MCs have been observed at distances other than 
at 1 AU, for example by the Helios spacecraft (e.g., Bothmer 
and Schwenn, 1998) and the Voyager spacecraft (e.g., Skoug 
et al., 2000; Burlaga et al., 2001). MCs have been revealed 
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in general to be large magnetic flux ropes in the solar wind, 
i.e., plasma embedded relatively strong magnetic fields of ap- 
proximately helical structure (e.g., Marubashi, 1986, 1997: 
Priest, 1990: Gosling, 1990; Lepping et al„ 1990; Burlaga, 
1995; Kumar and Rust, 1996). 

An interplanetary' coronal mass ejection (ICME) is an in- 
terplanetary manifestation of a coronal mass ejection (CME), 
which is strictly a solar event and well described by Schwenn 
(1996) and Gosling (1997); also see a comprehensive study 
of ICMEs by Richardson and Cane (2004) and Gopalsvvamy 
(2006). ICMEs are usually identified by their differences 
from the solar wind by some subset of about twelve po- 
tentially applicable quantities, including the same quanti- 
ties used for MCs (e.g., enhanced magnetic field intensity, 
smoothly changing field direction, and relatively low pro- 
ton temperature (e.g.. Gosling et al ., 1973) and low proton 
plasma beta), as well as by bidirectional streaming of elec- 
trons (e.g., Bame et al., 1981), bidirectional streaming of 
low energy protons (e.g., Marsden et al., 1987), high charge 
states of ions and compositional signatures, e.g., examination 
of Mgl0+/O6+ (e.g., Henke et al.,1998, 2001), low charge 
states (e.g., Zwickl et al., 1982; Galvin 1997; Burlaga et al., 
1998; Gopalswamy et al., 1998), singly charged He + (e.g., 
Schwenn et al., 1980), and bidirectional particle flows at cos- 
mic ray energies (~1 MeV) (e.g., Richardson et al., 2000), 
and bidirectional solar wind electron heat flux events (BDEs) 
(e.g.. Gosling et al., 1990), and including ground based data, 
the occurrence of a one or two step Forbush decrease (For- 
bush, 1937; Cane, 2000; Cane and Lario, 2006). Also see 
Schwenn (1996) who points out the important fact that of 
the many potential signatures of interplanetary' plasma clouds 
(meaning ICMEs, in our terminology) usually only a sub- 
set actually occurs for any given case. And see Russell and 
Shinde (2005) and Russell et al. (2005), who developed a de- 
rived parameter to define ICMEs which considers the com- 
bined plasma and field pressure perpendicular to the local 
magnetic field; this method strives to make ICME identifica- 
tion more unique. 

Both kinds of interplanetary structures (MCs and ICMEs), 
and others generally listed under the rubric of “plasma 
clouds,” where all of which often move faster than the nor- 
mal solar wind speed in their vicinity (at least at 1 AU), have 
been referred to as solar ejecta or just ejecta (e.g., Burlaga, 
1995). Ejecta are important because: (1) they always contain 
useful information about their dramatic birth conditions at 
the Sun (e.g., Webb et al., 2001; Gopalswamy et al., 2008), 
(2) they often drive interplanetary shocks which accelerate 
charged particles (e.g., Lee and Fisk, 1982; Kahler, 2001), 
and (3) portions of their extent, even their upstream sheath 
regions, are often markedly geoeffective (e.g., Burlaga et al., 
1981; McComas et al., 1988; McPherron, 1995; Wu et al., 
2006), because of generally having long-lasting strong south- 
ward magnetic fields and the usually higher than background 
solar wind speeds. Also see Luhmann (1997) who discusses 
the importance of CMEs to Space Weather generally. 


This study attempts to address the issue of why there are 
not more MCs observed within ICMEs as observed at 1 AU 
(e.g., for an early example of a MC associated with a coro- 
nal mass ejection see Burlaga et al., 1982). For example, 
it has been asserted that only 10-30% of MCs are observed 
within ICMEs (Gosling, 1990). Recent studies by Richard- 
son and Cane (2010) and Wu and Lepping (2010) agree with 
the upper limit of this range, although there is not complete 
agreement on this percentage in the space science commu- 
nity at large, e.g. see Xu et al. (2009); also see Richardson 
and Cane (2004). We suggest that the percentage of ICMEs 
containing MCs should be higher, perhaps much higher, be- 
cause an ICME is an interplanetary manifestation of a CME, 
and such an event is typically the result of a large explosive 
event on the Sun of twisted magnetic field lines (usually from 
an erupting prominence) and so also is a MC (Bothmer and 
Schwenn, 1994; Webb etal., 2000; Gopalswamy, 2006) . And 
the explosion itself is due to a magnetic process, although the 
details of which are not fully agreed on. Consistently, MCs 
depend on strong interplanetary magnetic fields for an es- 
sential part of their identification and so also do ICMEs, or 
they usually do. Hence, it is surprising that the occurrence 
of these two types of events are not more highly correlated. 
There may be several explanations. First, we point out that 
the definition of a MC being quite strict and that for a ICME 
being much less so is obviously part of the discrepancy in the 
occurrence rate of these two entities, and we speculate that if 
any magnetic structure even remotely resembling a MC were 
accepted as a real MC (e.g., accepting the greater number of 
magnetic cloud-like (MCL) objects found by Lepping et al., 
2005, in WIND data), then the number of MCs would obvi- 
ously increase dramatically. This may be an important, and 
perhaps the main, explanation for the discrepancy, especially 
since ICMEs correlate better with the combination of MCs 
and MCLs than with just bone fide MCs (see, e.g., Wu and 
Lepping, 2010). But below we suggest other explanations. 

When there is this "failure” to identify an existing MC 
within an ICME, is it at least occasionally due to a large 
closest approach (CA) passage of the spacecraft from the 
MC’s axis, making the MC more difficult to identify, i.e., is 
it due to a selection effect? Or is it a natural consequence of 
CME birth conditions, in the sense that apparently some (or 
many?) CMEs are simply not always born with the neces- 
sary smoothly twisted and strong magnetic field lines needed 
to indicate the existence of a MC? The possibility of the lat- 
ter explanation on a MC’s birth conditions will require study 
beyond our present scope, although we believe that it must 
be contributory for some cases. However, we may be able 
to shed some light on aspects of the geometrical argument, 
which will depend on actual and estimated |CA|s, and that 
will be our main focus here. (Note that we use the notation 
|CA| to indicate that we are not usually concerned with the 
sign of CA (= Y(j/Ro), which can be positive or negative, 
since Yo , the “impact parameter,” can be positive or negative 
depending on which “side” of the MC the spacecraft passes; 
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Ro is the MC’s estimated radius. Also CA can mean simply 
“closest approach,” as was used earlier, or relative closest ap- 
proach distance, which is \Yo\IRo- The context will usually 
make the usage clear.) 

It appears that this “geometrical” failure (i.e., one due to 
large |CA|) of MC identification is bound to occur, at least 
occasionally, especially for those cases where the spacecraft 
passes through the outer reaches of a MC where there is lit- 
tle change in field direction or intensity (If?!) and especially 
when \B\ is not too different from the ambient interplane- 
tary magnetic field. In those cases the MC would hardly be 
identifiable whether it is contained within an ICME or not. 
In this study we will, in fact, show quantitatively by using a 
large number of WIND MCs that there is a significant change 
in the variation of average field direction and average field 
intensity across a MC as the spacecraft’s |CA| grows from 
zero to a large fraction of the MC’s radius in a predictable 
way and becomes insignificantly small in both cases as |CA| 
approaches the MC’s radius, as might be expected. When 
there is a small change especially in field direction, along the 
spacecraft’s path within a MC, the MC is generally not going 
to be easily identifiable, and not likely amenable to analy- 
sis. For some of the above reasons, it is apparent that much 
of our analysis will require having relatively accurate esti- 
mates of |CA|, which is a quantity that is estimated by a MC 
parameter fitting model and whose accuracy is expected to 
depend on the sophistication of the model. Part of our anal- 
ysis, therefore, must show consistency among the estimated 
quantities of |CA|, Ro (also used as diameter, 2 Ro), Bo 
(the axial field intensity), and axial alignment, i.e., the rele- 
vant modeled quantities. Specifically, the field intensity part 
of the study will require the necessity to define a normalized 
magnetic field magnitude, and this is the reason for the use of 
the modeled parameter, Bo ■ And we must be able to find a 
proper model-dependent MC-coordinate system in which to 
do the analysis, the reason for knowledge of the axial align- 
ment, as well as be able to estimate sufficiently accurately 
|CA| itself. We will attempt to use the estimated quantities 
from the Lepping et al. (1990) MC parameter fitting model 
for a large number of cases from the WIND mission. How- 
ever, the study by Lepping et al. (2003) revealed that usually 
the poorest estimated quantity in the Lepping et al. (1990) 
model’s set of seven unknowns is CA, which plays a key 
role in this study. Hence, it is important that we take great 
care to prove its adequacy for the study. Further, we attempt 
to show the efficacy of using the estimated quantities Bo , 
Ro , and |CA| for their true values by calculating reasonably 
high correlation coefficients for the comparisons of expected 
normalized magnetic field intensity and directional changes 
vs. |CA| through so called “consistency” arguments, as ex- 
plained below. 

Our specific goals using WIND magnetic field data are to 
examine (and interpret) the following: (1) how severely non- 
uniform is the distribution of |CA|s for a large set of MCs 
(nonuniformity using a smaller number of cases was shown 
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Fig. I. Histograms of |CA| (= \ Yq\IRq \ (in %)), for various quality 
(Qo) subsets from WIND data, over the period from early 1995 to 
April 2006, confirming that there is a relatively severe drop off of 
identified MCs as a function of |CA| distance. 
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preliminarily by Lepping et al., 2006), and we will argue 
that this severity should not be expected, unless there is a 
selection effect mitigating against finding cases with large 
|CA|s, (2) parameter “consistency arguments,” described in 
detail below (concerning the relationships of Bo, Ro, and 
vs. |CA|), indeed do hold, providing useful and trustworthy 
estimates of these key parameters, especially |CA|, used crit- 
ically in this study, and finally (3) whether it is likely that 
actual MC sizes statistically control the sizes of estimated 
|CA|s, where all three of these questions are based on a large 
sample size of MCs. However we do not claim to know what 
the actual distribution of |CA| should be. 


2 The distribution of MC occurrences vs. |CA| 

We utilize the MC data from the WIND mission over the 
years from early 1995 to April 2006, and Table 2 on the Web- 
site http://lepmfi.gsfc.nasa.gov/mfi/mag_cloud_Sl.html pro- 
vides the model estimated values for Bo and jCAj (called 
\Y 0 IR 0 \ in %, in that table), based on the WIND magnetic 
field and plasma data (Lepping et al., 1995; Ogilvie et al„ 
1995). This table also provides Quality estimates (Qo= 1, 
2, or 3 for excellent, good, or poor, respectively, for each of 
these MCs); Appendix A of Lepping et al. (2006) gives the 
definition of Quality used here. The MCs will be split into 
various sets according to this Qo index. 

Here we show the distributions of |CA |s for various quality 
sets, given in Fig. 1 , i ,e., the distributions are based on the full 
set of WIND MCs, N - 98 and various subsets. (For the mo- 
ment we assume that the accuracies of the estimated |CA|s 
use here are high enough to generate a meaningful Fig. 1 . In 
Sects. 3.0 and 4.0 we defend this belief.) First, we point out 
that although [CA | (in %) should in principle go from 0 to 
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Fig. 2. The ideal cross-section of a magnetic cloud showing a space- 
craft passing through whose projected path is parallel to the Zq l 
axis in cloud (CL) coordinates (see text) with a closest approach 
distance Yq. The axis (normal to the plane of the figure) is 
parallel to the cloud’s axis and co-located. 


100%, we see that the distribution in Fig. 1 exceeds 100% 
for the <2o=3 set and hence also for the full set. This is sim- 
ply due to the fact that there is uncertainty in the estimations 
of the ICA|s for all sets (being especially poor for the Qo= 3 
set and even more so for the few cases where |CA| is near 
100%), as discussed by Lepping et al. (2003, 2004), and in 
the model there is no constraint that Yo must be equal to or 
less than Rq. But it is clear that in all cases we are consid- 
ering actual MCs, as has been discussed elsewhere (Lepping 
et al., 2006; and in this same paper a preliminary version 
of Fig. 1 is shown using a smaller number of MCs), so the 
actual |CA| must be less than 100%, except for special and 
rare cases discussed in the Summary and Conclusions where 
this overall issue is more extensively discussed. Therefore, 
the few cases where |CA| was erroneously estimated to be 
>100% by the Lepping et al. (1990) MC fitting routine, they 
are in actuality |CA| < 100% and usually are very near 100%, 
but we cannot give the exact actual value of |CA| , of course. 

We see in Fig. 1 that there is a somewhat dramatic drop 
of numbers of events as |CA| gets larger, for all subsets, be- 
coming for the N = 98 set, for example, only a small fraction 
of the total at [CA|=90%. A complete and accurate distribu- 
tion of actual MCs may not be exactly uniform, as we discuss 
in the Appendix A, but all of those in Fig. 1 seem to be unre- 
alistic and probably are due mainly to a selection effect as we 
get closer to a |CA| of 100%. In Fig. 1 we show each subset 
separately, because, even though the full set is apparently ad- 
equate for the argument, it has many poorly estimated |CAjs 
(mainly from the incorporated Qo = 3 set). But, even though 
the good sets, especially the Qo = 1 set, give more trustwor- 
thy results, there are far fewer cases (N- 23). In any case, 


all sets support the same basic conclusion: the distribution is 
markedly different from uniform. But the statistical accuracy 
of the estimates of the parameter |CA | used here is key to any 
conclusions drawn, of course. The following addresses that 
issue. 

3 Analysis: the ‘‘consistency arguments” 

Here we develop formulations for examining the variation 
in the magnetic field direction and intensity (|B|) vs. jCAj, 
separately, using the MCs with combined Quality of <2=1,2 
which provides 65 events. The main purpose here is to show 
that on average the model-estimated quantities Ro and |CA| 
are accurate enough for our analysis. There is usually little 
uncertainty in Ro . We will start by examining the change in 
observed field direction just w ithin the inbound and outbound 
boundaries for a set of MCs vs. |CA|. 

3.1 Angle a vs. |CA| 

We define an angle related to the maximum magnetic field di- 
rectional change across the MC for a given |CA|, after irrele- 
vant field “noise” is reduced by field averaging. For this pur- 
pose it is useful to define a so-called magnetic cloud (CL) co- 
ordinate system, where Xql is the unit vector along the MC’s 
axis, positive along the direction of the magnetic field at the 
axis, Y C l is in the MC’s cross-sectional plane along the clos- 
est approach vector, and Xcl x Ycl=Zcl; notice that Z cl 
is the projection of the spacecraft's trajectory in the MC’s 
cross-section, positive toward the Sun. The means of ob- 
taining the transformation matrix to go from geocentric solar 
ecliptic (GSE) coordinates (in which the magnetic field data 
originates in our case) to CL coordinates, as well as to see the 
final formulation of the matrix itself in terms of &a and (p,\ 
(which are the estimated latitude and longitude, respectively, 
of the MC’s axis in the GSE system) is given at the Website 
with URL http://Iepmfi.gsfc.nasa.gov/mfi/ediptic.html. 

The two angles 9a and <Pa uniquely define the CL coor- 
dinate system and have been estimated by the Lepping et 
al. (1990) MC fitting model for our analysis. Then we de- 
fine an angle a, which is ideally expected to be that between 
the magnetic field just inside the MC’s inbound boundary 
(Bin ) and.that just inside the outbound boundary (Bout), 
as projected into the MC’s cross-sectional plane (the Ycl- 
Zcl plane as seen in Fig. 2), by the following: 

a = cos -1 (Bin/ | Bin I x B OU t/|BoutI)cl- (1) 

We will use various interval averaging-lengths for these fields 
to help reduce field fluctuations that are considered “noise” 
in the sense that the fluctuations are not considered part of 
the underlying MC structure. The angle So is pictorially de- 
fined in Fig. 2. We claim that or ~ 2 So , where a is defined in 
terms of measured magnetic fields, from Eq. (1), and So is 
defined in terms of what is ideally expected of the magnetic 
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field just within the ideal MC's circular boundary. Notice 
that So is related to only the right side, whereas a depends 
on a difference between the right boundary ’s slope and the 
left boundary's slope, as is expressed in Eq. (1), and hence 
from this assumed symmetry we have the factor of 2. We 
then see from the expression for a circle, Y 2 +Z 2 =R 2 Q , and 
differentiation, that 

2Y o (dY/dZ)o+2Z o = 0, (2) 

where Eq. (2) is evaluated at point P, and hence, 
we use the subscript 0 on Yo and Zo ■ Since 

(dY/dZ)o=—tan8o=~sin8o/cos&Oi then Eq. (2) becomes, 

after substitution and dividing by —2Ro , the following 

Y 0 /Ro(sin&o/cos8o) = Z 0 /Ro- (3) 

We see from Fig. 2 that ideally sin 8o=ZolRo, and 
therefore, along with the definition of “Closest Approach” 
(|CA|) = \Yol\IRo, Eq. (3) becomes |CA|=cos/>o or, since 
a ^28o, it becomes 

|CAI cos(a/2), (4) 

where we note that ideally 0 < a <180°. (Notice that by as- 
suming symmetry we need only examine the right-hand side 
of Fig. 2, where for convenience we keep Zo > 0.) How well 
Eq. (4) is satisfied depends on how well the magnetic fields 
near the boundaries of the MC of interest satisfy the ideal 
circular model, i.e., one with low noise (after averaging the 
field) and no appreciable MC-expansion - and on how well 
Ro and |CA| have been estimated by the MC fitting model. 
Equation (4)’s validity will be checked using magnetic field 
data from actual MCs. 

3.2 Normalized magnitude of field vs. |CA| 

Likewise, we examine the variation in magnetic field inten- 
sity, | B | , but in the following way. For a given MC we use the 
estimated axial field from the Lepping et al. (1990) MC fit- 
ting model (called Bo), create an average of the actual field 
intensity (< \ B\ >) along the spacecraft’s path, and form the 
ratio < \B | >IBo- That is, we “normalize” the average field 
according to the axial field strength, where the average field 
value is obtained using the interval representing the space- 
craft’s path through the MC, i.e., the MC’s encounter’s dura- 
tion. In this way, we try to eliminate (or at least reduce) any 
problem with variation in “absolute” field strength among the 
set of MCs and obtain a ratio that is dimensionless and that 
will vary only with distance from the MC’s axis, i.e., with 
|CA|. Since there is a large variety of average field intensities 
across MCs, and similarly for the Bo s, we do not expect a 
good correlation for < jfi| > vs, (|CA|) orfor Bo vs.(|CA|), 
separately. This will be tested also. 

Now we ask, What do we expectfor < \B\ >/Bq vs. |CAj 
theoretically, if the Bessel function field of the constant a/ 
force free case for a MC used in Lepping et al. (1990) holds? 

www.ann-geophys .net/28/ 1 120 1 0/ 


We examine this question, noting that such a model is given, 
in CL coordinates, by: 

Ba = Bo Jo) tt-O (the Axial component, i.e., along Zcl 
in Fig. 2) 

Bj = BoHJi(w) (the Tangential component, i.e., 
perpendicular to A and R) 

Sr = 0, (the Radial component, i.e., along R in Fig. 2). (5a) 
where 

w — a'r , (5b) 

and where H is the handedness (H: +1 for right-handed and 
— 1 for left-handed); a ' is simply a proportionality constant 
that determines the pitch of the MC’s field lines (and not to 
be confused with the angle a used in the analysis above). 
Hence, the normalize (by Bo) magnitude of the field is 

\B\/Bo = f(u>), (6) 

where 

f(w) = J[(J 0 (w)) 2 + (Mw)) 2 ] (7) 

from Eq. (5). We now find the average of \B\/Bo (called 
[< 1 2? | > /fir; | theory) across the encounter-duration of the 
MC, based on this model. Here the average symbol < > 

represents an average across the duration, even though we 
will average over only the “first half” of the encounter and as- 
sume that the results hold for the full spacecraft’s encounter- 
duration, because of symmetry. So we integrate f{w) over 
dZ from 0.0 to Zo (where we drop the symbol “CL” on the 
coordinates as being understood) along the horizontal line 
of the spacecraft’s path as shown in Fig. 2, i.e., from the 
Z=0-point on the Y-axis to point P on the boundary (where 
Z-Zo)- Hence, the average of f(w) is 

[< |fi| > / B 0 ] theory — J f(w)dZ/ J dZ, 

(Z integrated from 0 to Zo; see Fig. 2), (8) 

We note that f(w) is a reasonably slowly varying function 
of w , and therefore it can be extremely well fit by a cubic in 
w: 

f(w)~ Ao + Aiw + A 2 W 2 + A 3 w 3 , (9) 

where A 0 = 0.99720, A^O.02412, A 2 =-0.17235, and 
A3=0.03296; the sigma for the cubic fit is 0.000865. Then 
we express [< |fi| >/fio]theory, which we call / for conve- 
nience here, as 

/ = £7/ (for J. 0 to 3), (10) 

where: 

lo - ( A 0 /Zo ) j dZ — Aq, and 

lj = (A j /Z 0 ) J w J dZ (for J: 1 to 3, 

and with integration limits of 0 to Zo), (11) 
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and where fdZ-Zo (in the denominator of Eq. 8), after 
f(w) from Eq. (9) is substituted into Eq. (8). We now re- 
late w and Z in the following way. First, we recognize that 

r 2 = Yl + Z 2 (12) 

along the spacecraft’s path (i.e., at fixed Yo), as shown 
in Fig. 2. Then, as done above, we note that specifically 
Y 2 o+Z 2 0 ^R 2 0 , or 

Zo = W( 1 - 1CA| 2 ), (13) 

(used for the upper limit in the integrals of Eqs. 1 1), and we 
notice that at the boundary' (where r=Ro) u>=2.40 (see Lep- 
ping et al., 1990), and therefore, a'=2.40 /Ro from Eq. (5b). 
Hence, using these facts and considering Eqs. (5b), and (12), 
we obtain 

w 2 = (2.40) 2 (|CA| 2 + Z 2 //?q) (14) 

Integrating Eqs. (11) using Eq. (14) yields (where again 
[< |5| >/f?o]theory= /): 

I = 0.99720 + 0.02895(1 + K) -0.33092(1 +2|CA| 2 ) 
+0.11390[l+3|CA| 2 (l + /O/2], (15) 

where 

AT^|CA| 2 [ln|(l+V(l-|CA| 2 ))/|CA||]/V(l-|CA| 2 ) (16) 

Note that in Eqs. (15) and (16) |CA| ranges over 0.0 to 
1.0; it is not in percentage. We will show below how well 
/ of Eq. (15) reproduces the actual WIND MC data for 
< \B\ >IBo vs. |CA|. And we stress that I was derived 
purely from the use of yj [(J o(w)) 2 +(J\(w)) 2 ] for the MC’s 
magnetic field magnitude, i.e., from the force free theory. If 
it does not agree with real MC data, then one or all of our 
estimates of |CA|, Bo, and the CL coordinate system (i.e., 
9a and 4>\ - see Sect. 3.1) would be suspect, as discussed 
above. 


4 Results 

In calculating < \B | > or < |R | >/Bo the averages are taken 
over the “duration” of each spacecraft passage for each MC 
from the basic field data in 1-min average form. But two 
types of averages, i.e., of interval-lengths of 15 min and 1 h, 
will be used for the computation of the angle a and cos(a/2) 
vs. |CA | ; this was done because we are focusing on sin- 
gle points just within the inbound and outbound boundaries 
for each MC, and such longer averages cause needed data- 
smoothing before computation. When all available MCs are 
used (i.e., for all levels of MC quality, Qo) we have yV=98 
cases, but this reduces to N= 65 cases when the poor cases 
(Qo- 3) are eliminated. 

In Table 1 we show what linear correlation coefficients 
(c.c.s) result when considering the angle a and cos(a/2) 
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Fig. 3. The angle a (black squares) and cos(a/2) (red *) plotted 
against jCA| (in %) all for Qo=l,2 (where iV= 65), for the 15 min 
average cases in the top panel (A) and for the 1 h average cases in 
the bottom panel (B). 


(from Eq. 1), both vs. |CA|, for different MC-quality con- 
siderations, in testing Eq. (4). In all cases as we go to greater 
quality (and therefore to a smaller number of cases), we ob- 
tain larger correlation coefficients. As we also see, the c.c.s 
for a and cos(a/2) are almost the same for the same Qo , but 
those for cos(a72) are consistently better, as expected, or at 
least the same (for the 1 h interval, N= 23). So in an aver- 
age way the assumption of the circular cross-section is quite 
good and supports “parameter consistency,” among a, Ro 
and |CAj. Figure 3 shows the results of a and cos(a/2) plot- 
ted against |CA| (in the same panel), all for Qo- 1,2 (where 
N =65) , for the 1 5 min average cases at the top (panel A) and 
for the 1 h average cases on the bottom (panel B). The 1 h 
set give the better c.c.’s for both a and cos(a/2), but slightly 
better for cos (a/2) (c.c. of 0.91), as expected. 

In Table 2 we show the c.c.s for < |J3| >, < |5| >/So,and 
So, all vs. | C A | . (Smax vs. |CA| was also considered with 
very poor results across the board. This is not surprising, 
since anomalous peaks in |S|, especially at the MC’s front 
are not uncommon. Hence, examining Smax is not useful in 
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Table 1. Linear correlation coefficients for a and cos(a/2) vs. |CA| . 


N of events 

Qo (and notes) 

15 min analysis interval 
c.c. {a} c.c.{cos(a/2)} 

1 h analysis interval 
c.c. {a} c.c.{cos(a/2)} 

98 

All events 

-0.71 

0.74 

-0.76 

0.80 

65 

1 and 2 

-0.82 

0.84 

-0.90 

0.91 

23 

1 only 

-0.93 

0.94 

-0.96 

0.96 


Table 2. Linear correlation coefficients for < |Bj >, < |B| >IBq and Bq vs. (CA|. 


N of events 

Qo (and notes) 

c.c.{< B >} 

c.c.{< B >IBo } 

c-c .{B 0 } 

98 

All events 

-0.27 

-0.93 

0.17 

65 

1 and 2 

-0.28 

-0.88 

-0.02 

23 

1 only 

-0.37 

-0.92 

0.00 


N= 65 Qg= 1 ,2 



ICA! 


Fig. 4. <\B\>, <\B\>/Bq, and Bq , in panels (A), (B), and (C), 
respectively, each plotted against |CA| and all for Qo = 1,2 (N= 65). 
The linear correlation coefficient (c.c.) for each quantity is shown 
in its panel. 


this type of analysis.) The c.c.s for < j B \ >iBo are all clus- 
tered near a good value of c.c .=—0.90 and, in fact, the one 
for “all events” (N=98) is unusually good considering that 
<2(9=3 cases are included. As we also expect, the c.c.s for Bo 
are either 0.00 or close to it. In fact, as we go to higher qual- 
ity the |c.c.| for Bo vs. |CA| gets lower, as it should, since 
there should be no correlation of Bo with |CA|. The c.c.s for 
< \B\ > vs. |CA| are poor, and this is also as expected. That 
is, without the normalization by Bo we get weak correlation 
of < \B\ > with |CA|. The good correlation that we see with 
<\B \ >/Bq is lost in examining only < \B \ >, no doubt be- 
cause there is a deleterious mixing of effects between actual 
dropping off of field intensity with |CA| and independent (of 
|CA |) variations of < \B\ > due to differing solar birth condi- 
tions and possible solar cycle changes. The main point here 
is that relative field intensity (< \B\ >/Bo ) behaves as ex- 
pected vs. |CA|, on average, strengthening the “consistency 
argument.” 

Figure 4 shows the results of < |B| >, < |B| >/Bo, and 
Bo plotted against |CA| in panels (a), (b), and (c), re- 
spectively, all al so for <2 o = 1 2 • We notice the ti ght cl ustering 
of < |B| >/ Bq with respect to ordinate-values in panel (b). 
As expected, the c.c.s for < \B\ > (panel a) and Bo (panel c) 
are very small and that for < |B| >lBo (panel b) is quite 
large and significant. We fit a quadratic curve to the data in 
panel (b) (now for [< |B| >/Z?o ]observations) and obtain the 
following: 

[< |B| > /B 0 ] obs' = Co + Ci |CAj + C 2 |CA| 2 

(for |CA|: 0.0 to 1 .0) (17) 

where C<9=0.8567, Ci=— 0.0304 and C 2 =— 0.3777, with a 
small sigma for the fit of 0.0389; the prime on 
[< \B\ >IBo lobs' indicates that the observations have been 
fitted to the quadratic curve. This can be thought of as sim- 
ply an overall quadratic smoothing of the data. In Fig. 5, 
panel (b) of Fig. 4 is again reproduced (i.e., [< |B| >IBq] 
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Fig. 5. Panel (B) of Fig. 4 showing [< \B\ >IBo ] (from ob- 
servations) vs. |CA| (in %), given by the black diamonds. The 
quadratic curve fitted version of this (== 7) is given by the dashed 
blue line (the primed quantity of Eq. 17). To this we superim- 
pose the [< |5| >/B0jtheory-curve (dotted red = /),from Eq. (15); 
see text. The solid red curve is given by Eq. (20); this is the 
[< |B| >/B(j]theory, mod estimate. The solid black curve at the bot- 
tom is given by (7 — /)//; it is a relative change in [< |B| >!Bq \ 
between observations and theory (before modification) - see Sect. 4. 

(from observations) vs. [CA| (in%), as black diamonds), as 
well as its quadratic fit (dashed blue curve, the primed quan- 
tity of Eq. 17) shown, and the [< |2?| >IBo] theory-curve 
(dotted red curve) , from Eq. ( 1 5) , is superimposed. It appears 
that, if Eq. (15) were multiplied by a constant ( 5 % 1.0) and 
slightly biased by constant “ a ,” we would nearly obtain the 
dashed blue curve in Fig. 5. Hence, 

[< IB] > /Bo]theory,mod — Sl + a. (18) 

How do we find 5 and a , and what do they mean? We now 
attack the fist question. In Fig. 6 we plot [< |B| >/Bo Itheory 
(= /) vs. [< |B| >IBo lobs' (= 7); i.e., each point is for the 
same |CA|. The result is nearly an exact straight line having 
a c.c. of 0.999. Hence, a linear fit of this line gives 

7 k =[< |B| > /Z?o] theory, mod = 5/ 4- a, (19) 

where we obtain 5=0.994 and a=+0.052 from the fit, with a 
small sigma of 0.0028. Hence, Eq. (18) becomes 

[< |£| > /B 0 ]theory,mod = 0.994/ +0.052, (20) 

which holds for all |CA|. The quantity 
[< \ B\ >/So]theory,mod from Eq. (20) is plotted in Fig. 5 as 
a solid red curve, and it very closely approximates the dashed 
blue curve in the figure which was a quadratic fit to the ob- 
served [< \B\ >/Bolobs ([< \b\ >!Bo Jobs') = 7)). But no- 
tice that to accomplish this good reproduction of observa- 
tions with theory vs. jCA|, we had to enhance the theoretical 
< jJS| >/Bo curve by a slight bias of +0.052 (after first mul- 
tiplying by 5=0.994). Even though “a” is not exactly zero, 
and 5 is not exactly 1 .0, we have achieved remarkably good 
agreement between theory and observations, using our model 


c.c, = 0.999 



[<IBI>/B„] 

THECFt 

Fig. 6 . [< \B\ >IBq fobs' (from the dashed blue curve of Fig. 5) vs. 
[< |B| >/Bq] theory (from the dotted red curve of Fig. 5). These 
two sets are very well correlated, with c.c.=0.999. (This good cor- 
relation enables us to calculate S and a in Eqs. 18 and 19). The 
diagonal dashed line has a slope of exactly 1 .0 and shows, by com- 
parison, how closely the data follows this slope. 


|CA|s, Bq s, and derived coordinate system (through the val- 
ues for 9 a and 4>a - see Sect. 3.1), suggesting again that we 
can have confidence in these estimated entities on average. 
But the differences between theory (before modification) and 
observations (as fitted by the quadratic) are worth investigat- 
ing quantitatively, so we now form the normalize difference 
ratio 


Ratios (7 -/)//, (21) 

This Ratio is a measure of the relative change in 
[< |B| >IB 0 ] between fitted-observations and theory vs. 
|CA|, and is given by the black curve in the bottom of Fig. 5. 
Notice that this Ratio is relatively low and steady, suggesting 
that our relative measure of the effects of solar wind inter- 
action and expansion on the structure of a MC (i.e., outside 
of the simple constant a' force free solution to a MC’s field 
structure) is, on average, weak and almost independent of the 
radial distance from the MC’s axis. However, there is a mod- 
est increase in Ratio as |CA| increases, possibly suggesting 
that there is a greater influence of external compression in 
the outer reaches of the MC than in the inner regions, as one 
might expect. This finding is an interesting by-product of this 
analysis, but since there are many sources of uncertainly in 
its derivation, we should be cautious in our interpretation. 
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N= 65 c.c.= 0.28 



CAI 


Fig. 7. A plot of the estimated diameters of the MCs (2 Ro) vs. |CA| 
(in %) for the combined Qo= 1,2 cases (N= 65), where c.c.=0.28, 
and where each quality-set is shown separately; for Qo=\ (N= 23. 
by * in blue) and for Qo=l (N= 42, by diamonds in red). 

5 |CA| vs. R 0 

In Fig. 7 we show a plot of the estimated diameter, 2 Ro, 
vs. |CA| for the WIND MCs of quality sets Qo = 1 (blue *) 
and Qo - 2 (red diamonds). The c.c. is 0.28 for the two sets 
combined; separately the c.c.’s are 0.11 (Qo= 1) and 0.38 
(Qo=2), both being quite small. See Table 3 which com- 
pares the c.c.’s for 2 Ro vs. |CA| for various combinations of 
quality sets - again all showing small c.c.’s. It is evident that 
there is no significant linear correlation between the |CA|s 
and the sizes of the MCs. So it is not likely that the MC sizes 
statistically control the sizes of |CA|s. It may be interesting 
to note, however, as also shown in Fig. 7, that for the Qo- 1 ,2 
set of MCs the largest and smallest events tend to have Qo - 2 
quality evaluations. 

6 Summary and conclusions 

WIND magnetic field data were used in this analysis cov- 
ering the period from early 1995 until mid April 2006, and 
providing N - 98 MCs. Results are categorized according to 
the Quality (Qo) of the events considered. Appendix A of 
Lepping etal. (2006) gives the definition of Qo - Most of the 
analysis was carried out using MCs of Qo=l and 2 combined 
(where Qo spans values of 1 (excellent), 2 (good), and 3 
(poor)), i.e., the poor cases were not used for the main results, 
but some of their consequences were shown. We explained 
that we occasionally show estimated values of |CA|s greater 
than 100%, even though rarely are they actually so. In prin- 
ciple no value of |CA| for any realistic MC should be greatlr 
than 100%, unless the MC is a core-annulus case (e.g., see 
Lepping et al„ 2006), and they are uncommon; in this case 
the (inner) core’s radius is at Ro • Even for the poorly fitted 


Table 3. Linear correlation coefficients for 2Rq vs. iCAJ . 


N of events 

<2 o ( an d notes) 

C.C. 

98 

All events 

0.19 

65 

i and 2 

0.28 

23 

1 only 

0.11 


MCs (all of which are Qo = 3 cases), we choose to present the 
model’s actual output estimate of |CA| rather than artificially 
assigning a value of |CA| ~ 100%, even though the actual 
value must be near 100%±(£), where the error E (i.e., E 
with one sigma uncertainty) is usually quite big for many of 
the Qo = 3 cases. That it is even possible for model-estimated 
values of |CA| to be greater than 100% is due to the fact that 
the model does not constrain Yo to be equal to or less than 
Ro (where recall that |CA| = \Yo\/Ro)', they are both free 
parameters in the model. In summarizing a large early set 
of WIND MCs, Lepping et al. (2006, Fig. 13) show the typ- 
ically large spreads of sigma{|CA|}s on MCs generally. In 
many cases, even when Qo is 1 or 2 (i.e., the good cases) 
the range of E ' s on |CA| is large, i.e., from 15% to 55%, and 
for the Qo= 3 cases the E's on |CA| can range from 20% 
to 70%, or even larger (but rarely). So, for example, if the 
estimated |CA| is 120±40%, this surely allows for a |CA| 
smaller than 100% (i.e., within the range of 80% to 100%), 
but we cannot give more information on where it is exactly in 
that range. In fact, we have previously shown that the |CA| 
parameter (i.e., Yo) is usually the poorest determined of the 
full set of seven parameters that we fit with our model for 
any quality classification; see Lepping et al. (2003, 2004) on 
a study of the fit parameters for the force free model. In any 
case, the MC examples with poor parameter fittings (<2o=3) 
are of minor concern to us here, because we rely mainly on 
the <2=1,2 cases in the quantitative aspects of this study. It 
is precisely because of our concern about the usefulness and 
accuracy of estimated |CA| by our model, and |CA|’s impor- 
tant role in this study, that we emphasized the consistency 
arguments of Sect. 3. 

Magnetic field observations within MCs indicate that as 
|CA| (in %) approaches 100%, both the magnitude and di- 
rection profiles of the field (vs. time) are typically somewhat 
flat, making it hard to identify a MC based only on magnetic 
field observations when there is a distant spacecraft passage. 
(However, examples of such flattening are not shown explic- 
itly in this work.) This is true, of course, whether the MC 
is “within” an 1CME or not. Also the success of our “con- 
sistency” arguments (Sects. 3 and 4 concerning both the Inl- 
and a -profiles) aids us in trusting sufficiently, and in an aver- 
age way, the model-estimated quantities that were important 
to the analysis. These quantities are Ro , |CA|, Bo , and the 
MC-coordinate system in which the analysis was carried out. 
(Again we stress that the “consistency” arguments hold best 
for the Qo= 1 ,2 cases and only for a large number of events 
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(N- 65). For the Qo= 3 (poor) cases we should not always 
expect such satisfying results, as the correlation coefficients 
in Table 1 show for the a study, for example.) That being 
determined, we generate the distribution of the number of 
MCs vs. |CA| and showed in Fig. 1 how it decreases as |CA| 
increases, whether the event is ICME-related or not, where 
again a large number of cases are considered and for various 
quality subsets. Figure 1 shows that there is a dramatic drop 
off of the numbers of events as |CA| becomes larger for all 
subsets; Fig. A1 shows similar results where the MCs were 
parameterized according to high vs. low inclinations, split at 
45° (shown also to be true if split into three parts, i.e., at 30° 
and 60°). In reality a complete distribution of |CA|s may not 
be exactly uniform, but the estimated |CA|s in Fig. 1 seem to 
have an unrealistically rapid drop-off, which we argue is due 
in part to a selection effect. By this we mean that it is usually 
more difficult to recognize a MC when the |CA| is large and 
where the magnitude of the field is relatively low and does 
not change very much in magnitude or direction over the du- 
ration of the spacecraft passage. In fact, it has been observed 
that | B | at large |CA| within a MC can drop to approximately 
ambient values especially in the rear of the MC (not shown 
in this work). 

We first examined the field direction change from the 
inbound- to outbound-boundary within an average MC, 
through the “a-analysis,” and showed a relatively good corre- 
spondence between cos a/2 vs. |CA| statistically, as expected 
ideally for our model (see Figs. 2 and 3). The correspondence 
was especially good when carried out using the 1-h averages 
of the field which yielded a c.c. of 0.91 , as shown in the bot- 
tom panel of Fig. 3. We then concentrated on field intensity 
and normalized intensity. 

In our attempt to theoretically reproduce (from the force 
free model) the actual average fall-off of the normalized 
magnetic field intensity in [< |B| >/fio]theory vs. |CA|, we 
had to slightly decrease the theory-curve for 
[< |B| >/Bq | theory by a factor of 0.994 (derived from S) 
and slightly bias it by a =+0.052; see Fig. 5 and the associ- 
ated text. Although the values of 5 and a represent small 
adjustments, they are not insignificant. We are not sure that 
we fully understand their cause. However, they are likely, 
at least in part, the result of the fact that MCs interact with 
the solar wind which usually causes field compression within 
the MC, especially in the outer reaches. (But to be complete, 
the effects of MC expansion, especially in explaining the dis- 
placed location of Bmax within a MC, for example, should 
be folded into this analysis also - see Farrugia et al. (1993). 
MC expansion occurs for most cases at 1 AU.) Field com- 
pression is often present in the upstream part of the MC 
but not so often in the downstream part. But the effect of 
the overall interaction (and expansion) is to cause an aver- 
age intensity increase over the full duration, especially in the 
ZcL-component of B . and therefore also in < \B\ >, without 
much change in the direction of the actual field from the sim- 
ple force free field direction. And both effects (expansion and 


compression) usually tend to cause a very noticeable asym- 
metry in the |B (-profile across the MC such that the field is 
usually stronger in its early part, and sometimes markedly so, 
at least at 1 AU. (But also at 1 AU some relatively small num- 
ber of MCs are quite different, i.e., some have relatively sym- 
metric |B|s, peaking near the center, some |B|s are nearly 
flat, and some even peak somewhat late in the MC.) It is in- 
teresting to note that at greater distances from the Sun most 
MCs are seen to grow even larger according to a simple ex- 
pansion rule, as seen in Voyager observations, for example 
(Burlaga and Behannon, 1982; Burlaga, 1995). However, 
even at distances of 2 to 4AU the field magnitude within 
a MC often retains some asymmetry in \B\, but to a lesser 
extent than at 1 AU (Burlaga and Behannon, 1982). Even- 
tually, at much greater distances, the expansion is expected 
to significantly decrease or disappear and with it most of the 
asymmetry in | B\ as the MC’s magnetic field merges into the 
background interplanetary magnetic field. However, the de- 
tails of this are not completely clear, so it is an area that needs 
further investigation 

“Ratio” ( Y - 1)1 1 (see Eq. 21) is a measure of the relative 
change in [< |J?| >IBq ] between observations and theory as 
a function of |CA[, and its examination is another way of 
looking at the impact of S not being exactly 1 .0 and “a” not 
being exactly 0.0. Our analysis shows that Ratio is small 
and relatively steady with respect to small or medium sized 
|CA|s (i.e., up to ~50%), but it slowly increases as |CA| ap- 
proaches 100%, as shown by the black curve in Fig. 5 (bot- 
tom); actually something close to a plateau is reached. We 
interpret this to mean that the interaction/expansion effects 
on the structure of a MC are on average somewhat weak and 
almost independent of the radial distance from the MC’s axis 
but are more noticeable in the outer reaches, due to external 
effects. Ratio’s behavior is an interesting by-product of this 
study. (This behavior says nothing about (< B | >/Bo)’s re- 
lationship to longitude within the MC, of course.) We stress 
that these results do not apply for any specific MC, where 
interaction/expansion effects may be strong and very notice- 
able, or very weak, depending on specific conditions. 

Regardless of the well-known shortcomings arising from 
the simplicity of the Lepping et al. (1990) MC fitting model 
(force free of constant a'), we find that key fit-parameters are 
internally consistent on average especially for the < 20 = 1,2 
WIND set. And this set of MCs represents about 2/3 of all 
cases know for that mission (i.e., up to and including the year 
2006). Those key parameters are Ro, Bo, and CA, as well 
as the latitude ( 9a ) and longitude (</>&) of the MC’s axis that 
provide the magnetic cloud coordinate system in which we 
carried out the analysis. And since handedness, another fit 
parameter in the model, is always known, we are accounting 
for the relatively good accuracy on average of 6 of the 7 fit 
parameters for the <2o=l ,2 WIND set. This is another by- 
product of our analysis, and probably of greater importance 
than the one mentioned above. 
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We suggest that a spacecraft passing through MC at a large 
(CAj can lead to the MC-encounter being overlooked, be- 
cause MCs are much harder to identify when the |CA| is 
large. And this may happen frequently. However, this may 
not completely answer our original question (on the occur- 
rences of ICMEs vs. MCs), and we must look for other fac- 
tors, such as peculiarities of CME birth conditions. And 
a thorough analysis of the subject would require MC-solar 
source correlations. But such correlations are difficult and 
sometimes made all the more so, because some originally 
expelled MCs (seen at the Sun by any means) are expected 
to be missed by Earth (e.g., see Gopalswamy et al., 2008). 

This point should be stressed: ICMEs are defined in terms 
of a far larger number of potential physical indicators than 
are used for MCs (e.g., see Schwenn, 1996), meaning per- 
haps 12 vs. 4 parameters, respectively, i.e., a ratio of 3 times 
as many for ICMEs than MCs. But different subsets of these 
12 physical parameters are utilized for different ICMEs mak- 
ing it easier not to miss identifying an ICME encounter. And 
we point out that most of the ICME identifiers are indepen- 
dent of |CA|, or at least weakly dependent on it, in contrast 
to the situation for MCs. This may be the principal cause of 
a poor MC-ICME correlation (e.g.. Gosling, 1990, and Wu 
and Lepping, 2010), but specifically for the following rea- 
son. Suppose that a MC/CME structure is expelled from the 
Sun such that, as it passes a spacecraft at Earth, with say the 
encounter being far off the axis of the internal enhanced \B\ 
portion of the overall helical field structure (indicating a poor 
or no MC encounter), but nevertheless with some of the 12 
ICME-parameters being measured indicating the presence of 
an ICME. Many such ICME identifiers do not depend on a 
parameter such as |CA| (e.g., bi-directional streaming elec- 
trons), and more than a few ICME-indicators may be present. 
Hence, the MC-ICME correlation will be driven down for 
such cases. But, if the probability of observing both MCs 
and ICMEs depended on average |CA| in the same way, re- 
gardless of the particular value of (CAj for any case, corre- 
lation for a large set would be quite high, provided that the 
structures were highly correlated at the Sun in the first place, 
of course, which we do assume. This is a key element of our 
argument as to why, at least partly, an unusually low num- 
ber of MCs are identified within ICMEs as observed at 1 AU 
(^30% or smaller); see Gosling (1990) and Wu and Lep- 
ping (2010). And we speculate that if any magnetic structure 
even remotely possessing a MC field geometry and intensity 
were accepted as a real MC in the identification process, then 
the ratio of MCs-to-ICMEs would probably increase dramat- 
ically. 

Russell and Shinde (2005) and Russell et al. (2005) dis- 
cuss the well known problem of non-uniqueness associated 
with the identifications of ICMEs. They developed a new pa- 
rameter to define ICMEs whereby the net plasma and field 
pressure perpendicular to the local magnetic field is used, 
with consideration of its difference from the ambient value, 
and even in their work several classes of events appear to 


Table Al. Coefficients for Eq. (Al) a . 


Set considered 

N 

K 0 

b 

a 

All Qo=l2 

65 

47.9 

1.8xl0“ 6 

40.7 

0°<!#aI <45° 

43 

47.0 

-8.6xl0” 8 

40.9 

45°<!e A l <90° 

22 

49.1 

8.1 xlO- 7 

40.7 


a for jCAj measured in %. 

be possible. The Russell et al. (2005) approach to identify- 
ing ICMEs is perhaps the most satisfying yet in an attempt 
to avoid ambiguities. But, in any case, the broadness of the 
ICME’s definition compared to that for MCs is very likely 
not the only factor in accounting for the lack of apparent co- 
incidence of these two types of structures. 

Finally, we see no significant correlation between the es- 
timated ]CA| values and the sizes of the MCs (measured by 
the ideally estimated diameter, 2 Ro by our model) for the 
WIND’S Qo=l2 MC set; see Fig. 7. This was shown to be 
true also for the full set of N = 98 MCs, but displayed here for 
only the Qo= 1,2 set (N=65). So, in general, it is unlikely 
that the distribution of sizes of MCs statistically controls the 
distribution of the observed sizes of |CA|s. 

Appendix A 

What should the distribution of |CA| really look like gener- 
ally, if there were no selection effect? By using only obser- 
vational data from a single spacecraft (or even from multiple 
spacecraft, if the number of them or their positions are in- 
sufficient) we may never know the answer to this question. 
However, it may be revealing to examine again the distribu- 
tion of |CA|s for the WIND set of MCs but now parameter- 
ized according to whether or not each MC is very' inclined 
or not, representing a combination of birth conditions, birth 
locations, and resulting solar wind interactions. Hence, we 
choose to look at those cases of 1 6 a I (absolute value of the es- 
timated latitude of the axis of a MC) between 0° and 45 c and 
those between 45° and 90°, separately, where 9 a is strictly 
defined in Sect. 3.1. We show the results in Fig. Al. As in 
Fig. 1 (showing the observed distribution of |CA[) these his- 
tograms show that there is a severe drop off of identified MCs 
as a function of |CA| for both sets, and they look similar. We 
try' a least-squares fit to each of the three distributions (where 
the center point of each bar is the value fitted) to the normal 
distribution of variable amplitude, i.e., to 

N(x) = K 0 G + bx)exp(-x 2 /2a 2 ) (Al) 

where x = |CA| (in %). We considered Eq. (Al) to be sym- 
metric about CA=0 and allowed both positive and negative 
CA for fitting purposes, providing 8 points to fit for each 
curve, with the understanding that only the positive CA do- 
main is to be retained. Kg,b,a and the number ( N ) of MCs 


www.ann-geophys .net/28/ 1 / 20 1 0/ 


Ann. Geophys.,28, 1-14,2010 



12 


R. P. Lepping and C.-C. Wu: Selection effects in identifying magnetic clouds 



ICAI 


Fig. Al. Histograms of |CA| (= \Yq\IRo\ 0 n %)). for Qo = 1 an< J 
2 combined split into those cases of \6&\ between 0° and 45° and 
those between 45° and 90° separately, as well as for the full set of 
N = 65 cases. Curves of normal distributions are least-squares fitted 
to all three distributions, showing very similar properties. 

in each set are given in Table Al . Notice that b is extremely 
small in each case and can be ignored, so the simple normal 
distribution is adequate for fitting these three distributions. 
We then see that all curves are essentially the same (within 
small errors) giving Ko ~ 48 and cx ~ 41 . We carried out 
this same exercise with a split of three intervals of the |0 aIs 
according to 0°-to-30°, 30°-to-60°, and 60°-to-90° and ob- 
tained essentially the same results as for the 0°-to-45° and 
45°-to-90° split. So it does not matter if we split the distri- 
bution into low and high inclination MCs; we will get the 
same distribution in either case. We believe that this tells us 
that the specific inclinations of the MCs, regardless of their 
source, probably do not contribute to any discrepancy be- 
tween the distributions of |CA[ for MCs vs. ICMEs. 

Supplementary material related to this article 
is available online at: 
http://www.ann-geophys.nct/28/l/2010/ 
angeo-28- l-2010-supplement.pdf. 
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